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Model for reflection and transmission matrices of nanowire end facets 
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Nanowires show a large potential for various electrooptical devices, such as light emitting diodes, solar cells 
and nanowire lasers. We present a method developed to calculate the modal reflection and transmission 
matrix at the end facets of a waveguide of arbitrary cross-section, resulting in a generalized version of the 
Fresnel equations. The reflection can be conveniently computed using Fast Fourier Transforms, once the 
waveguide modes are known. We demonstrate that the reflection coefficient is qualitatively described by 
two main parameters, the modal field confinement and the average Fresnel reflection of the plane waves 
constituting the waveguide mode. 
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I. INTRODUCTION 

Nanowires have received significant attention due 
to promising applications in electrooptical devices. 
Nanowire lasers^, light emitting diodes^ and solar 
cells^ have already been demonstrated. Semiconductor 
nanowires are characterized by large aspect ratios and 
high index contrast. Various techniques can be used to 
accurately control the composition and geometry of the 
nanowires. Some work has been done to analyze the opti- 
cal properties of semiconductor nanowires^"—, including 
modal properties, modal gain and optical emission. More 
knowledge is however needed to better understand and 
predict the optical response of the nanowires, in order 
to improve the nanowire devices. For nanowire lasers, 
two especially important properties are the reflection and 
transmission at the end facets. The reflection is impor- 
tant for the laser resonator, whereas the transmission 
yields the far field radiation. 

The reflection properties of nanowires have been in- 
vestigated by Maslov et alA and Henneghien et al^, us- 
ing a 3d finite element method, and also by Hugonin et 
ali^ using coupled- wave analysisi^. With this work they 
have contributed with important knowledge concerning 
the reflection properties of nanowires. Calculation of 
the reflection and transmission at an interface between 
a waveguide and some half-space is however a 2d prob- 
lem connecting the true waveguide modes at each side of 
the interface; thus it is not necessary to perform a full 3d 
analysis. In addition, the cross coupling between differ- 
ent modes, including coupling to radiation modes, and 
the transmission far field, are not obtained directly with 
the finite clement methods. 

In planar dielectric waveguides, various techniques 
have been used to analyze step discontinuitieai^"— . Sev- 
eral of the techniques are based on matching modal fields 
at the boundarjii^r— . Most attention has been given to 
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conventional waveguide structures with low index con- 
trast. 

In this work we use the electromagnetic boundary con- 
ditions to obtain matrix Fresnel equations. These equa- 
tions directly give the modal refiection matrix [ry], and 
the transmission matrix [tim]- The model has mainly 
been intended for calculations on nanowires, but is also 
useful for other waveguide geometries. A key aspect of 
nanowires are their small lateral size, often on the scale 
of the light wavelength, leading to significant diffraction 
effects. In our model, these diffraction effects are mani- 
fested as cross coupling to other modes, including radia- 
tion modes. The model is more compact and simple in its 
form than the previous Id techniques, and it is directly 
applicable to waveguides of arbitrary cross-section, both 
in Id (planar waveguides) and 2d. 

This article is organized as follows. The general re- 
fiection model for an arbitrary geometry is discussed in 
Section II. In Section III the special case of a planar 
waveguide is presented. The simplifications involved in 
the planar case facilitate the illustration of some impor- 
tant aspects of the general model. Results for hexagonal 
nanowires are presented in Section IV: a wide bandgap 
nanowire based on a material similar to ZnO or GaN, a 
narrow bandgap pure GaAs nanowire and hetcrostruc- 
tured GaAs nanowires with Alo.3Gao.7As claddings. Fi- 
nally, in Section V, the model is used to calculate the 
refiection properties of a multimode optical fiber. 



II. CALCULATING THE REFLECTION AND 
TRANSMISSION MATRICES 

We will now demonstrate how the electromagnetic 
boundary conditions may be reformulated as matrix Fres- 
nel equations. The waveguide is oriented along the z-axis 
and terminated at z = 0, see Fig. [TJ No particular geom- 
etry is assumed for the cross-section. At each side of the 
boundary, z = 0, the mode spectrum is discretized. The 
discretization can be done using e.g. periodic boundary 
conditions, i.e. the single waveguide is approximated by a 
periodic array of nanowires. We surround the waveguide 
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tions, (l2cl) and ()2bp follow from ([2a|) and (j2dp . respec- 
tively. The modes are orthogonal in the sense 
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Here, the area of integration is defined by —L < x < L 
and —L < y < L. The unit vector in the z-direction is z, 
and 
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FIG. 1. Waveguide oriented along the z-axis with end facet 
at 2 = 0. The waveguide core and cladding have refractive 
indices rico and rici, respectively. The ambient medium for 
2 > is assumed to be homogeneous with refractive index ria. 
The fields are discretized using periodic boundary conditions, 
i.e. the structure is assumed to be periodic with periodicity 
2L in both x and y direction, where 2L ^ 2p . 



with a quadratic box of size (2L)^, and let the area within 
this box constitute the unit cell of the periodic grating. 
The period in both x- and y-direction is thus 2L. The 
parameter L should be large enough, so that neighboring 
waveguides do not perturb the field of the bound waveg- 
uide modes. In the ambient half-space z > there is a 
homogeneous isotropic medium with refractive index na. 
The electric field of the modes can be written 



Smix, y) = £m exp^ikxX + ikyy). 



(1) 



The magnetic field, Hm, is described in the same way. 
Each mode is labeled by a collective mode index m = 
fcy, pol). Specifying the wavevector and the polar- 
ization. The polarization, pol, is either TE or TM. The 
real transverse wavevectors are = s-nd ky = qj^, 
where p and q are positive integers. The modal propa- 
gation constant k^ is given by k^ = n^w^/c^ — k^ — ky, 
where is the refractive index of the half-space z > 0, 
and c is the vacuum light velocity. The four constant 
vectors can be expressedSS 
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A 



•-m(TM) 



ykxkz , kykz,- 



(2a) 

{kl + kl)), (2b) 
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where A=\ j y (fc^ -I- A:^)|fcz|2L2 ^ £^ is the permittivity, 
and /i is the permeability. From Maxwell's curl cqua- 



For the waveguide [z < 0), the modes can be di- 
vided into two sets: Discrete bound modes, and con- 
tinuous radiation modes^i. Due to periodic boundary 
conditions (or alternatively, metallic boundary condi- 
tions), however, the continuous set of radiation modes 
is converted into a discrete set. We will accordingly de- 
scribe the modal fields by the discrete set {e^, hi}, where 
Gi = ei{x, y) and — hi(x, y) are the electric and mag- 
netic fields, respectively, of mode i. 

The transversal electric field must be continuous at 
the boundary. Assuming incoming mode ei{x,y) we can 
write 



(t) _ 



, Hm'-'m ^ 



(5) 



valid for all a; and y. Here is the reflection coefficient 
from mode i to mode j, and tim is the transmission co- 
efficient from mode i (in the waveguide) to mode m (in 
the half-space z > 0). The superscript (t) stands for the 
transversal component (x and y components) of the vec- 
tor. The boundary condition for the transversal magnetic 
field is similarly 
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Eqs. (O and ^ can be combined as follows. Take the 
vector product between (O and 'H^*{x, y), and integrate 
over the unit cell. Similarly, take the vector product 
between £^*(x,y) and ([6]), and integrate over the unit 
cell. This yields 
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Eliminating 



we obtain 



(9) 



Eqs. and dH) are multimode, matrix versions of 
the electromagnetic boundary conditions. The matrices 
[^im] and have first index i (designating the waveg- 
uide mode) and second index m (designating the mode 
of the homogeneous medium). Solving ([9]) yields the ma- 
trix form of Fresnel's equation to calculate the reflection 
matrix: 



[nj] = ([$™] - [*»,„])([$»,„] + [^un]y 



(10) 



In (jlOp . we have assumed that there are at least as many 
m values as the number of i's; thus the equation can 
be solved with the help of Moore-Penrose pseudoinverse 
(denoted by superscript +). Once the reflection matrix 
is found, the transmission matrix can be calculated with 
the help of ©. 

The matrices and VPim can be written out in terms 
of the Cartesian components of the modal fields as fol- 
lows: 
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Here we have defined 
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cell 



The other field components are defined similarly. 

The distance between neighboring transverse wavevec- 
tors must be sufficiently small, so that a superposition of 
the associated modes in the half-space z > can describe 
the modes of the waveguide (z < 0) accurately. In other 
words, we require the spacing 



Akx < -, Aky < -, 



(13) 



where 2C is the maximum extension of the bound modes. 
Incq. (HSl) will be fulfilled for L > £. 

For practical computation, we first need to compute 
the waveguide modes {e^, h;}. Once the modes are avail- 
able, the inner products ([8]) are found, using ([TT|) - ([T2l) . 
The inner products can be conveniently computed using 
the fast Fourier transform. 

In pT7)) we are free to multiply each side by any factor 
dependent on m — m(pol, fc^., fey). Ideally this will not 
affect the solution. However, in practice the result after 



pseudoinversion will differ somewhat depending on this 
normalization, as the different equations in the system 
are weighted differently. 

Given a solution [rij] it is interesting to calculate the 
error, measured as the deviation from perfect match of 
the tangential fields. We assume incoming mode i, and 
denote the deviation of the transversal electric and mag- 
netic fields by Ae^*' and A.hf\ respectively. We can 
define corresponding transformed quantities as follows: 



Ai?„„ = f^i / Aef)xHW*(^,y).2dA (14) 
|K(m)| 2 J 

cell 

(Ai?™)* ^ / £^\x^y) X Ahf> .zd^. (15) 

|k(to)| 2 J 

coll 

Denote the transmission coefficients given by ([7a|) and 

(E9 by t^j2 a-iid i-^^ respectively. Taking [t^J^l the 
definition of the transmission matrix, the electric field is 
perfectly matched. The error in the tangential magnetic 
field is then given by the deviation between the left and 
right hand side of (|7b)) : 



|Ai/„„l - 



(e) _ ^(Ti) 



(16) 



Similarly we find the error of the electric field when 
is taken as the definition of the transmission matrix: 



lAK;, 



(e) _ Ah) 

1 r>i ^ ■» 'yft 



It is therefore natural to define the error 



(17) 



(18) 



In this work we have chosen to work with a homoge- 
neous ambient medium, due to the simplification that the 
inner product can be written as Fourier transforms when 
the modes {£m, T~Lm\ are complex exponentials. Eqs. ([7])- 
are however generally valid, and can also be applied 
when the transmission medium is inhomogeneous, as long 
as one is able to find the modes of the medium. Further- 
more, Eqs. ©-(O are directly applicable if one chooses 
to use metallic boundary to discretize the field on the 
waveguide side, i.e., let a metallic box of size (2L)^ sur- 
round the waveguide. 

If the dimension of the waveguide is significantly 
smaller than the material wavelength of the light, the 
waveguide will be highly diffractive. In this case, the 
waveguide facet can be viewed as a point source, spread- 
ing spherical wavefronts. The spherical wavefronts that 
propagate along the plane z = 0, are reflected at the 
artificial boundary, and may thus couple back into the 
waveguide. To prevent such artificial refiections from dis- 
turbing the solution, we may introduce some loss, i.e., 
the permittivity everywhere is transformed according to 
1769. (Here Eq is the vacuum permittivity.) The 
loss parameter 7 should be small enough not to alter the 
reflection properties of the boundary significantly. An 
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indication of this is obtained by comparing the Fresnel 
reflection coefficients for reflection at an interface with in- 
dex contrast nco/n^, to those of an interface with contrast 
V '^co + *7/ + *7- The Fresnel coefficients should re- 
main approximately unaltered for all angles of incidence. 

The modal fields and eigenvalues, of the waveg- 
uide arc given from the boundary conditions and the 
Helmholtz equation applied to e.g., the longitudi- 
nal field component . 



e!^) = 



(19) 



Here n{x,y) = rico in the core region and n(x,y) = 
in the cladding region. Note that by introduc- 
ing the same loss in the core as in the cladding re- 
gion (n{x,y)'^ — n{x,y)'^ +*7), the Helmholtz equation 
(|19P remains unaltered if we perform the transformation 



III. PLANAR WAVEGUIDE 



2a 




FIG. 2. Planar waveguide with thickness 2a oriented along 
z-axis with end facet at z = 0. The waveguide core and 
cladding has refractive indices tIco and rid, respectively. The 
ambient medium for z > is assumed to be homogeneous 
with refractive index n,. 



We will here consider the special case of a planar 
waveguide. This special case is interesting because the 
simple geometry enables us to find both radiation and 
bound modes analytically, and reduction from 2d to Id 
offers better probabilities to test convergence criteria. It 
also has interesting applications, being the geometry of 
conventional semiconductor laser diodes. 

A planar, dielectric waveguide with thickness 2a is lo- 
cated along the z-axis, see Fig. [21 The refractive index 
of the core (|a;| < a) and cladding > a) of the waveg- 
uide is rico and rid. The waveguide terminates at z = 0; 
for z > there is a homogeneous nonmagnetic medium 
with refractive index ria. We divide the modal set into 
two sets depending on the polarization (TE and TM). 
In this waveguide any TE (TM) polarized mode can be 
described using a superposition of TE (TM) plane waves 
([2]). According to this yields 'i'im, = ^im = when 
i and m describe orthogonal polarization. Thus the two 
sets of modes will not couple to each other through the 
boundary z = 0; so we may consider them separately. 
Here we will only consider TE modes. 

To discretize the modes on each side of the boundary, 
we impose periodic boundary conditions: The fields are 
supposed to be periodic along the a:-axis; for example, 
the fields at x = +L correspond to the fields at x = —L. 
The artificial period 2L satisfies 2L ^ a. 

The electric field of the modes for z > can be written 



CP 

"--TE 



A^y^uj^ cxp{ik^x + ifcf z)y. 



(20) 



Here = pn/L, (fcf )^ — {ns.uj/cf — (fcg)^, and p is an 
integer labeling the mode. The transversal component of 
the magnetic field is 



niP 
''-TE 



APkP 



cxp{ikPx + ifc^z)x. 



(21) 



The normalization constants 
such that 



l/^y\kz\L arc chosen 



idx 



kP* 



Opp' . 



(22) 



For the waveguide (z < 0), the modes are a little bit 
more complicated. Due to reflection symmetry about the 
plane x = 0, the set of modes can everywhere be divided 
into even modes and odd modes. These two sets will not 
couple to each other, so we can consider them separately. 
We take the z-dependence to be exp(i/3z), and solve the 
wave equation in each region < a; < a and a < x < L. 
By requiring an even mode, the electric field in the re- 
gion < X < ais B cos{k'^x), where k'^ = w^n^^/c^ — 
and i? is a real normalization constant. For a < x < L 
the electric field can be written C sm{kxX + ip), where 
k^ = ^^ci^^/c^ ~ and C and ip are real constants. The 
electromagnetic boundary conditions require the electric 
field and its derivative with respect to x to be continuous 
at a; = a and x — L. Moreover, the periodic boundary 
condition combined with even symmetry requires that 
dsin{kxX + fi)/dx\x=L — 0. Taken together, these condi- 
tions require that 



(23) 



k'^ tan(A;^a) — —k^ tan[fc^(L — a)]. 



Denoting the j-th solution of the dispersion relation ((23)) 
by superscript j, we can write the modal fields as follows: 



-il}j{x) exp(i/3^z)y. 



(24) 



where 



\ cos{k'Jx), 
cos[ki{L - x) 



for 0<x<a, ^25) 
for a < a; < L, 
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and 



cos[ki{L — a)] 



(26) 



The transversal component of the magnetic field is given 

by 



(27) 



The normalization constants are chosen to obtain 
orthonormality, similarly to 



1 



ej X h* ■ zdx 



1 133* 



ip*{x)'iljjx)dx 
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W\ 



^33' 



(28) 
(29) 



confirm our numerical results. Secondly, the core refrac- 
tive index is quite similar to the refractive index of GaAs, 
which will be considered in the next section. 

The necessary loss parameter, 7, can be reduced (in- 
creased) by increasing (decreasing) the half-width L ac- 
cordingly. Note however that increasing L quickly be- 
comes numerically challenging for a 2d structure, as the 
density of modes scales with 1/L in Id and in 2d. 

For comparison, we have used two sets of simulation pa- 
rameters; both L = 60 • a, 7 = 0.1 and L = 240 • a, 
7 = 0.025. The latter set of parameters describes the 
real structure more accurate, whereas the first set of pa- 
rameters have significant advantages in terms of compu- 
tational resources. Comparing the Fresnel reflection coef- 
ficient at normal incidence for an interface with refractive 
index contrast 7ico/«a with and without loss, the differ- 
ence is |Ar|/|r| « 0.032 for 7 = 0.1, and |Ar|/|r| « 0.008 
for 7 = 0.025. The waveguide was described using all 
modes with real B. in addition to the evanescent modes 



The matrices and [$,p] can be calculated using ® : ^^^^ ^""^^ > -(10c^a/c)2. For half-width L = 60, this 




ijji exp(— ifc^a;)dx 



(30a) 
(30b) 



From (jlOp . the matrix Fresnel equation for a planar 
waveguide is 



(31) 



Note that for plane waves, Eq. ([31]) is reduced to the stan- 
dard Fresnel reflection coefficient for TE polarization, 
''TE = {k\ — k1)/{k]. + k1), where k\ is the z-componcnt 
of the propagation constant in medium i. 

The transmission coefficients are found from ([7]): 



corresponds to 577 modes at auj/c=3, and less for lower 
frequencies. With this high number of modes, we ob- 
tain a very good field match at the boundary. However, 
accurate results for the reflection matrix are obtained 
with fewer modes. For example, increasing the cut-off 
to — (5ti;a/c)^, i.e., reducing the number of modes by ap- 
proximately a factor of 2, the reflection matrix is roughly 
the same, while the field match at the boundary is worse. 

Fig.|3]shows the reflection of the fundamental TE mode 
as a function of ua/c. The figure also shows ri2, i.e. the 
fraction of the fundamental mode reflected into the flrst 
excited mode. The result corresponds well with results 
from Vahidpour et ali^. In the geometric optics limit, 
the modes can be described as plane waves with incident 
angles (f>^ = arccos(/3^/(ncow/c)). The reflection can thus 
be described using the Fresnel equations. This limit is 
indicated in Fig. [31 i.e., the Fresnel reflection coefficient 
of a plane wave with index contrast rico/na, and angle 



Up = I f^p + ^T-y^-jp 



(32a) 



(32b) 



Similarly to the general case in Section II, we may 
introduce artificial loss by modifying the refractive in- 
dices according to n1 — )• -I- 17, n^^ — >■ n^^ -|- 27, and 
n^j — >■ n^j -I- 17. For the homogeneous medium, this im- 
plies a modification (fcf)^ (^f)^ +*7(^/c)^, and for 
the waveguide {(3^)^ — (/3^)^ + ij {uj/cf . The transverse 
wavenumbers thus remain unaltered. 

We now consider a planar waveguide of width 2a and 
refractive index n^a = -s/TO, completely surrounded by 
vacuum. The motivation for this choice of rico is twofold. 
Firstly, Vahidpour et alr^ have previously performed cal- 
culations for this structure, giving us the possibility to 



Without the use of the artificial absorption, we ob- 
served oscillations in the reflection as a function of fre- 
quency. We attribute these oscillations to artificial reflec- 
tion from the hard-wall boundary conditions. Different 
frequencies will give rise to varying waveguide- mirror res- 
onances. By adding the artificial absorption, these oscil- 
lations were removed. The good agreement with Vahid- 
pour et alr^ verifies that the absorption has not altered 
the reflection properties significantly, even for 7 — 0.1. 
The agreement between the calculations with 7 = 0.1 
and 7 = 0.025 is also very good. In the following, we 
will look closer into the results from the calculation us- 
ing 7 = 0.1. 

Fig. m shows the transmission coefficient of the funda- 
mental mode at wa/c=3, calculated as the mean of the 
transmission matrices obtained from the two equations 
(|32ap and (|32bp . The deviation (fTH]) between the trans- 
mission resulting from the two equations is also shown. 
The deviation increases at the spatial frequency corre- 
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FIG. 3. The reflection coeflicient as a function of auj/c for 
a planar waveguide of width 2a and refractive index Wco ~ 
y/W. The figure shows the fraction of the fundamental mode 
amplitude reflected back into itself, |rii|, and the fraction 
of the fundamental mode amplitude reflected into the first 
excited mode, |ri2|. Two sets of parameters have been used, 
one with 7 = jh = 0.1, and one with j = jl ~ 0.025. 



spending to the highest order waveguide mode included 
in the model. Fig. [5] shows the match of the electric and 
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FIG. 4. The transmission coefficient from the fundamental 
mode into plane waves with wavevectors k^. The mismatch 
(II8II is also shown. 



magnetic transverse fields e*^*^ and h^*\ with the funda- 
mental mode incident auj/c=i. The match is very good, 
especially for the dominating spatial frequencies, i.e. the 
kx values where the transmission is the largest. The 
match can be improved further by increasing the number 
of modes in the calculation, i.e. by decreasing the cut-off 
value of Re/3^. 

It is interesting to compare the modal reflection co- 



FIG. 5. The transverse fields e'*' and /i'*' at the boundary 
2 = 0, the mismatch between the fields at z = 0"*" and z = 
0~ is also indicated. Z = ^/JloJs^ is the impedance of the 
ambient medium. 



efficients with the Fresnel coefficients. From Fig. [31 we 
see that the fundamental mode was well described by the 
Fresnel equations for aw/c > 2. To investigate how the 
behavior is for higher order modes, we calculate the re- 
flection matrix at auj/c = 10. The same modal cut-off 
as previously was used, the loss was 7 ~ 0.1, and the 
half-width of the unit cell was reduced io L = 10a. The 
increased normalized frequency, auj/c, leads to reduced 
diffraction effects, thus enabling a smaller L parameter. 
Fig. [6] shows the diagonal elements of the reflection ma- 
trix, as a function of the modal angle (j)^ . 
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FIG. 6. The refiection coefficient as a function of the modal 
angle (j)-' = arccos(/3-'/(ncoij.'/c)). For comparison the Fresnel 
refiection coefficient and the confinement factor rj are also 
shown. 

Two main factors affect the rcffcction in the diffraction 
regime. Firstly, the field is only partly confined within 
the waveguide, and the field outside the waveguide region 
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experience a reduced refractive index contrast of ricxln^ 
rather than nco/jT-a- This effect will be more pronounced 
with decreasing confinement, i.e., for higher order modes. 
The confinement factor 77, defined as the fraction of power 
within the core, is included in Fig. [6] for reference. Sec- 
ondly, each mode can be described as a set of plane waves, 
centered around the angle ^ . Decreasing aw / c yields an 
increasing angular spread of each mode. The effect of 
this angular spread is a smoothing of the reflection spec- 
trum. In Fig. [ni it is seen that around the critical angle 
for total internal reflection the reflection is smeared out 
compared to the Fresnel reflection. Furthermore, for de- 
creasing confinement factor the reflection is seen to de- 
crease compared to the Fresnel reflection. 

IV. HEXAGONAL NANOWIRES 

We will now focus our attention towards two kinds 
of semiconductor nanowires, narrow band gap nanowires 
based on ZnO or GaN, and wide band gap nanowires 
based on GaAs. These structures have been reported 
to have a hexagonal cross-section. Investigations have 
previously been done to compare the reflection properties 
of hexagonal and circular nanowires^; they were reported 
to be approximately the same when comparing structures 
with the same cross-section area. It is thus convenient to 
describe hexagonal nanowires using an effective radius, 
Peff^, defined so that a circle with radius peff has the same 
area as a hexagon of sides a, i.e., a « 1.1 x peff- In our 
calculations we will use the exact hexagonal geometry; 
Poff is however chosen to describe the nanowire. This 
facilitates comparisons with previous work^ and circular 
nanowires. 

The modes of the hexagonal nanowires were found us- 
ing the commercial software, Comsol Multiphysics"^'^ . 
For these hexagonal structures it was more convenient 
to use a metallic rather than a periodic boundary condi- 
tion to discretize the waveguide modes. 

The extension to a 2d cross-section as well as the 
strong diffraction present in the nanowires, impose some 
additional challenges compared to the planar case. To 
maintain a feasible model in terms of computer power 
and time consumption, some care was taken to limit the 
number of modes needed to describe the field accurately. 
Firstly, evanescent modes were ignored. Secondly, as the 
wavelength decreases the diffraction becomes less pro- 
nounced and the radial extension of the field will de- 
crease. The half-width of the unit cell, L, was thus set to 
be frequency dependent, decreasing linearly as a function 
of wavelength according to Llo/c = constant. 

A. Wide band gap nanowires 

We consider a hexagonal nanowire with refractive in- 
dex TT-co = 2.45, surrounded by vacuum. This value 
is close to refractive index of ZnO (2.33)^ and GaN 



(2.60)2^ at room temperature. Such a nanowire has been 
studied previously by Maslov et alA and Henneghien et 
al^. ZnO and GaN have wide bandgaps, and will typi- 
cally lase around 380 and 375 nm2^, respectively, at 
room temperature. Note that ZnO and GaN are slightly 
anisotropic, but this has not been taken into account. 

The simulation parameter for the artificial absorption 
was set to 7 = 0.1, and the half- width of the imposed 
unit cell was L = 20pcff at the lowest frequency, and 
decreasing as a function of increasing frequency. 

Fig. [7] shows the diagonal elements of the reflection 
matrix for a nanowire with rico = 2.45, for the 3 modes 
with the lowest cut-off. The next excited mode has its 
cut-off value at peSOJ / c ~ 1.5; a higher number of modes 
should be included in the calculation in order to describe 
the reflection properties of such higher order modes accu- 
rately. We used around 220 waveguide modes to describe 
the fleld in the nanowire. The result coincide quite well 
with results reported earlieri^, especially for the funda- 
mental mode (HEii). The reflection coefficient deviates 
slightly more for the TEqi and TMqi mode, especially 
for the TMoi in the frequency regime where this mode 
has a very low confinement. 

After a regime where the reflection is dominated by 
the reduced modal conflnemcnt, each mode approach the 
Fresnel reflection as observed in the planar case. The 
fundamental mode is in this region close the Fresnel re- 
flection at normal incidence for a refractive index contrast 
2.45/1, i.e 0.42. For the higher order modes TEoi and 
TMqi the average incident angle of the plane waves con- 
stituting the modes will be higher, so that in the geomet- 
ric optics picture one would have total internal reflection. 
The reflection is however limited by the reduced conflne- 
mcnt. This effect explains the fact that the reflection of 
these modes rise quickly, to exceed the reflection of the 
fundamental mode. Highest possible reflection coefficient 
and single mode operation are important to optimize a 
laser. This is achieved for a radius around 60 nm for 
the nanowire considered here, i.e., just below cut-off of 
the first excited mode. In this region the modal reflec- 
tion coefficient of the fundamental mode is 0.35. For 
larger radii, the TEoi mode will start to dominate the 
modal spectrum due to its much higher reflection coeffi- 
cient, whereas for smaller radii the refiection coefficient 
will decrease due to limited confinement. 

If one is interested only in the qualitative behavior of 
the fundamental mode, it is possible to obtain a quick 
and quite good estimate of the reflection coefficient us- 
ing only bound modes. The result of such an estimate is 
shown in Fig. |S1 The resemblance to the more accurate 
calculations (Fig. [7|) is quite remarkable for the funda- 
mental mode. Note that all that is needed is to Fourier 
transform less than 10 modes. One might simplify further 
by approximating the nanowire as a cylindric waveguide; 
in this case the modes could be described using analytic 
expressions from standard textbooks^!. 
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FIG. 7. The reflection coeflicient for the lowest order modes 
for a hexagonal nanowire with rico = 2.45 surrounded by vac- 
uum, rici — 1. The superscript designate the degenerate fun- 
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FIG. 8. Rough estimate of the reflection coefficient for the 
lowest order modes for a ZnO hexagonal nanowire surrounded 
by vacuum, using only bound modes. 



B. GaAs-based hexagonal nanowires 

GaAs has a more narrow band-gap compared to ZnO, 
and GaAs will typically lase around 850-870 nm at room 
temperature. As will be seen, nanowire lasers will thus 
need to have a larger diameter compared to the wider 
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FIG. 9. Reflection coefficient of the three lowest order modes 
for a series of GaAs based nanowires, with core refractive in- 
dex rico = 3.3. The cladding consists of Alo.3Gao.7As shells, 
rici = 3.15, of varying thickness td- The waveguides are sur- 
rounded by vacuum. 



band gap nanowire considered in the previous section. 
Surface defects arc often a problem for nanowires due to 
the high surface to volume ratio of nanowires. To reduce 
the amount of surface defects, it might be advantageous 
to grow a radial shell surrounding the nanowire. Such 
a shell will serve as a cladding material of the nanowire 
waveguide. A radial shell of AlGaAs^^ or GaAsP— is 
commonly used. We have investigated how the reflection 
properties vary with increasing cladding thickness for an 
Alo.3Gao.7As shell. Fig. [5] shows the reflection for one 
pure GaAs nanowire, and three GaAs/Alo.sGao.yAs core- 
shell nanowires, all surrounded by vacuum. We have used 
refractive index values for infrared wavelengths, i.e. the 
core refractive index is rico ~ 3.3, while the refractive 
index of the cladding is rid = 3.15. 

The artificial loss parameter was again taken to be 
7 = 0.1. For the pure GaAs nanowire, the factor Pes/L 
was 20 for the lowest frequency, linearly decreasing for 
increasing frequency. The nanowires with cladding had 
the same cell sizes as the pure wire. To better enable 
comparison between the two nanowire systems, we have 
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used a scaled abscissa, ricoPcffW / c, taking into account the 
difference in refractive index. Comparing Fig. [7] and tlie 
pure nanowire in Fig. [9l we see tliat the quahtative be- 
havior is very similar, but that the higher refractive index 
of GaAs leads to an increased reflection. In the geomet- 
ric optics limit, the fundamental mode of the pure GaAs 
nanowire should approach the Frcsncl reflection at nor- 
mal incidence, which is 0.53. From Fig. [HI we see that the 
fundamental mode reflection is close to this value from 
ncow/cpeff = 2.5. To optimize the lasing performance, 
one should operate the laser at a wavelength where the 
nanowire waveguide is single mode, and with the highest 
possible fundamental mode reflection. It is also desir- 
able to achieve lasing at the smallest possible diameter. 
Operating the nanowire laser at 850 nm, the effective 
radius of a pure GaAs nanowire should be around 100 
nm. As could be expected, the increase in reflection com- 
pared to ZnO or GaN based nanowires is proportional to 
the Fresnel reflection at normal incidence, 0.53/0.42, and 
the optimal radius has increased with a factor equal to 
the ratio of material wavelengths at the lasing frequency. 
As the cladding thickness increases we see a clear trend 
in terms of reflection properties. Higher reflection can 
be obtained for smaller radii, but the reflection curve 
of the fundamental mode also has a higher gradient in 
the single mode region. It thus becomes more difficult 
to design an efficient nanowire laser, due to shrinking of 
the single mode region where the reflection of the funda- 
mental mode is high. Another feature that is important 
when choosing the cladding thickness is the fraction of 
the modal field confined in the active GaAs cor o^'^° , to 
maximize gain. 

We now consider the pure GaAs nanowire in the region 
around the single mode cut-off. If some of the design pa- 
rameters deviate slightly, one might excite several modes. 
It is therefore useful to know how the modes couple to 
each other in this region, and how excitation of higher 
order modes can be diagnosed by analyzing the far field. 

Fig. [TU] shows the reflection matrix at Ticow/cpcff = 2.5, 
i.e., just above the single mode cut-off. The reflection 
matrix figure shows the power reflection for the bound 
modes, [l^ij ; the last element of each row/column dis- 
plays the total reflected energy from/to the given bound 
mode to/from radiation modes. From the reflection ma- 
trix we see that the coupling between the various bound 
modes is quite weak, whereas there is a significant cou- 
pling to radiation modes. Cross coupling between bound 
modes would have strengthened the competing excited 
modes at the expense of the fundamental modes, and 
thus e.g. deteriorate the beam profile. 



C. Far field 

The beam profile in the far field is an important as- 
pect of the lasing performance. The far field radiation 
of a circular ZnO nanowire laser has previously been in- 
vestigated numerically by Maslov et al^. An advantage 
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FIG. 10. The power reflection matrix [Ir'ijP] for the bound 
modes of a GaAs hexagonal nanowire, with peff = 60 nm 
and Ucoi^/cpefi = 2.5. The last element of each row/column 
displays the total reflected energy from/to the given bound 
mode to/from radiation modes. 

of our method is that the far field is directly obtainable 
from the transmission matrix, so no further simulations 
are necessary. 

For the calculations of the reflection and transmis- 
sion, as mentioned previously, we utilized an artificial 
absorption to avoid refiections from the imposed bound- 
ary conditions. The reflection and transmission coeffi- 
cients of the actual structure were thus approximated by 
the reflection and transmission coefficients of the struc- 
ture with loss 7. Once [r^j] and [tim] have been found, 
one should again consider the actual structure. The ar- 
tificial absorption should thus not be included in the far 
field calculations, i.e., the medium for z > is isotropic 
and homogeneous with real refractive index ria. The 
far field is given from the Fourier transfornt^^ of the 
field at z = 0; which can be conveniently found from 
the transmission matrix [tim]- For incoming mode i, 
the Fourier transform of the electric field at the inter- 
face is simply i]i{ka:,ky) = £,n{k^,ky,TE)ti,ni{k^,ky,TE} + 

£ni{k^,ky,TM)ti,m.{k^,ky,TM), and similarly for the magnetic 
field. As a measure of the beam profile we calculate the 
z-component of Poynting's vector. In the far field, at the 
point specified by direction [k^, ky,kz], this power flux is 

\Sf{kx,ky)\ CX kl (|ii,m(fc^,fc„.TE)P + l^i,m(fc^,fc„.TM)P) • 

(33) 

We investigated the far field of the pure GaAs nanowire 
at ricoPoffw/c = 2.75, using the simulation data from 
Fig. [HI Due to periodic boundary conditions, the far 
field becomes very coarse. We therefore performed an 
interpolation by padding with zeros in the Fourier do- 
main to restore sufficient spatial resolution. Fig. llll shows 
\S'^{kx,ky)\, when the fundamental mode is incident to 
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the end facet of the pure GaAs nanowire. Any presence 
of higher order incoming modes, in addition to the fun- 
damental mode, win lead to a change in the far field. 
For reference, the far fields for incoming modes TEqi 
and TMoi are therefore shown in Fig. [121 The qualita- 
tive behavior is in accordance with previously reported 
data^. Note that upon reflection of i.e. the fundamental 
mode, higher order modes will also be excited and give 
contributions to the far field. Fig. [11] is thus not simply 
the far field of the fundamental mode, but of the total 
sum of modes excited upon reflection of the fundamen- 
tal mode. Also note that the hexagonal cross-section is 
visible in Fig. [11] If the fundamental mode is excited 
in the nanowire, the far field beam profile will be Gaus- 
sian like, with a peak intensity in the forward direction 
(kx = ky = 0). The two higher order modes will on the 
other hand lead to zero emission in the forward direction. 
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FIG. 12. Far fields, \S?^\, for the pure GaAs nanowire excited 
by the modes TEoi and TMoi; ricoPcffi^/c — 2.75. 



have here used L = lOp, and no artificial absorption was 
introduced. Furthermore, we have neglected coupling to 
radiation modes; this may cause a reduced accuracy for 
the highest excited modes. The modes of the optical fiber 
were found using the analytical expressions for cylindrical 
waveguides^!. The diagonal elements of the reflection 
matrix for light at the wavelength 633 nm are shown in 
Figs. [T^KHl The reflection is plotted as a function of 
the modal angle (/)■' = arccos(/3^7(ncow/c)). 



FIG. 11. Far field, 15*^1, for the pure GaAs nanowire excited 
by the degenerate pair of fundamental modes; ricoPeH'jj/c = 
2.75. 



V. MULTIMODE OPTICAL FIBER 

Reflection properties of multimode optical fibers are 
interesting in order to understand the modal composition 
of the light after reflection. We have calculated the facet 
reflection of two commercially available silicon step index 
optical fibers from Thorlabs, with rico = 1.4570 and rici — 
1.4537 at the wavelength 633 nm. The core radii of the 
two fibers are p = 12.5 /im and p = 5 pm. The cladding 
thickness of both fibers is 125 ^m, being so thick that we 
have approximated the claddings to be infinite. These 
multimode fibers have Lup/c ~ 125 and ujp/c 50, so 
the diffraction effects are less pronounced than those of 
the structures we have studied previously. Therefore the 
unit cell can be reduced compared to the nanowires. We 



0.189 



0.187 



0.185 



0.183 



I 0.181 



0.179 



0.177 



0.175 







□ 
* 



TM 
TE 

Hybrid 
Fr.-TM 

Fr.-TE 



/ * 



0.05 



0.1 



0.15 



^[rad] 



FIG. 13. The diagonal reflection matrix elements, \rii\, for 
the bound modes of a silicon step- index optical fiber of radius 
p = 12.5 pm, terminated in air. 
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waveguides, and more specifically nanowire waveguides, 
where the end facet reflection is vital for the nanowire 
lasing performance. The end facet reflection was studied 
for ZnO nanowires, and for GaAs/Alo.sGao.rAs nanowire 
heterostructures with various cladding thicknesses. In 
general the reflection of the fundamental mode is signifi- 
cantly different from the plane wave Fresnel reflection for 
normal incidence. However, for both systems the max- 
imum reflection in the single mode region was found to 
be close to this value. As a by-product, the model pro- 
vide the field at the boundary, so that the far field beam 
profile can be directly calculated. The model is also well 
suited for calculations on multimode fibers, as we have 
demonstrated for two commercial fibers. For all waveg- 
uide geometries, we demonstrate that the reflection from 
a waveguide end facet is mainly governed by two factors, 
the modal field confinement and the Fresnel reflection of 
the plane waves constituting the waveguide mode. 



FIG. 14. The diagonal reflection matrix elements, |rii|, for 
the bound modes of a silicon step- index optical fiber of radius 
p = 5 /im, terminated in air. 



The Fresnel reflection coefficients tte and ttm for 
plane waves with incident angle <\>' , are also shown in 
Figs. [TBlfni Note that in analogy to the observations 
for the planar waveguide, the reflection coefficient of the 
modes with TE, TM polarization resemble the tte, ''tm 
coefficients, respectively. The polarization properties of 
the hybrid modes can be seen as a mixture of TE and 
TM polarization, and they have reflection coefficients in- 
termediate between the tte and rxM coefficients. The 
comparison to the Fresnel reflection is best for the lowest 
order modes (with smallest angles 0^). Similarly to the 
planar case this can be explained by the reduced con- 
finement of the higher order modes. Note however that 
the accuracy is somewhat decreased for the highest order 
modes, due to the neglect of radiation modes. 

For these optical fibers, the intermodal coupling is 
ver} - ! weak. Thu maAimiuii coupling hum the fmidamcii- 
tal mode into the other modes, max [ / [ ri^ l ), for 
these two optical fibers was 1.3 x 10"'^ and 4.1 x 10"'* 
for the fiber with p = 5 /im and p = 12.5 /im, respec- 
tively. As expected, there is more modal cross coupling 
in the smaller fiber, where diffraction effects are more 
pronounced. 



VI. CONCLUSION 

We have presented a general formalism for calculating 
the end facet reflection in various waveguide geometries. 
The method is a natural generalization of the Fresnel 
equations to corresponding matrix equations, describing 
the reflection and transmission into all modes. A spe- 
cial focus in this article has been on highly diffractive 
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